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It is described how the standard Poisson bracket formulas should be modified in order to incorporate integrals 
of divergences into the Hamiltonian formalism and why this is necessary. Examples from Einstein gravity and 
Yang-Mills gauge field theory are given. 



1. INTRODUCTION 

Hamiltonian mechanics traditionally serves as 
an etalon part of the mathematical physics both 
for physicists and mathematicians ||l|J^. Many 
of its geometrical constructions are exported now 
to the field theory. For example, the Schouten- 
Nijenhuis bracket |^ turned out to be extremely 
useful in the search for integrable models during 
the last 20 years. In the pioneering article [|| 
this approach had been called as the formal varia- 
tional calculus because in proving theorems it was 
possible to replace functionals by functions, even 
polynomials. Presently there are some mono- 
graphs available where the method is presented 
comprehensively 

The main restriction required by this approach 
is the freedom to integrate by parts. This require- 
ment is fulfilled if the fields are rapidly decaying 
at spatial infinity (as the massive fields usually 
do) or if the periodic boundary conditions are 
imposed. However, these cases do not cover all 
types of the physically interesting boundary be- 
haviour. Often one has to deal with a slow decay 
of the gauge or gravitational fields at spatial infin- 
ity 1^,^ , with nonperiodical boundary conditions 
for a media moving in a finite domain |^ or for 
integrable models pO|-p^]. Can we still use the 
Hamiltonian approach in such cases? Of course, 
the answer is positive and the cited papers just 
do that. What we are trying to do here is to pro- 
pose a more general framework for the nontrivial 
boundary problems. 

In this report we intend to show that the main 
geometrical concepts of the formal variational cal- 



culus, and, consequently, the Hamiltonian me- 
chanics, can be preserved without the standard 
requirement on the freedom to integrate by parts. 
In other words we propose an approach that does 
not neglect surface integrals. This purpose can be 
achieved by introduction of a new grading into the 
formal variational calculus |13| and a new pairing 
compatible with this grading. First of all, the re- 
vision consists in the modification of formulas for 
the Poisson bracket in field theory p^ . 

2. THE POISSON BRACKET 

Briefly then, we will modify the standard for- 
mula for the Poisson bracket by surface terms in 
such a way that all its general axiomatic prop- 
erties (bilinearity, antisymmetry, Jacobi identity 
and closeness, i.e., the requirement to remain in 
a given space of functionals) ||^ will be preserved 
without discarding any term in the course of in- 
tegration by parts. 

Being one of the main elements of the Hamil- 
tonian formalism the Poisson bracket itself is not 
elementary and may be considered as a composite 
structure 



{F,G} = dGjdi^J*. 



(1) 



Its elements are: the differentials of functionals 
dF, dG, the Poisson bivector 5* and the pairing 
(or interior product) operation _|. It turns out 
that in order to take care of all surface integrals 
we should revise all the three constituents of the 
bracket. 

There are two ways to write a local functional: 



as an integral of some smooth function 
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of the fields and their spatial derivatives up to 
some finite order over the prescribed domain 
in R", or as the integral over all the space K" but 
with the characteristic function of the domain 9n 
included into the integrand 
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elusion 
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In turn this leads us to the con- 
that transformations of the form 



(9) 



(2) 



Henceforth we will consider the space M", use the 
Einstein rule for summations and the multi-index 
notations J = (ji, j„) where ji > 



\J\^jl + -+Jn- (3) 



The reader unfamiliar with these notations may 
first have in mind the one-dimensional case, then 
J simply is the order of spatial coordinate deriva- 
tive. Binomial coefficients for multi-indices are 



Jn 
kn J 



1 otherwise. 



(4) 



(5) 



are, generally speaking, canonical only up to sur- 
face terms. Finally, as the standard proof of 
the Jacobi identity in mechanics is based on the 
commutativity of the mixed second derivatives, 
in field theory the Jacobi identity for functionals 
may not be true even if the fields themselves have 
the canonical Poisson brackets. 

To improve the situation we allow arbitrary 
variations on the boundary 



dF = 



o<PA Jn 
Dj iEi{f)6<p, 



f'AiHA) 



(10) 



where the differential is written consequently by 
using the full variational derivative 



0(pA 



(11) 



With the help of them we introduce the so-called 
higher Eulerian operators p5 16 ,^ 



where 



(.J+i) 



D{\..Di- 



(6) 



(7) 



In the framework of the standard approach the 
differential of a local functional is given by the 
Euler-Lagrange derivative E\{f) 



the Frechet derivative 
df 



D, 



(12) 



and the higher Eulerian operators (|^). 

The second constituent, a Poisson bivector, is 
given, loosely speaking, by Poisson brackets of the 
fields. These brackets are called local if they are 
proportional to the (5-function and a finite num- 
ber of its derivatives (ultralocal, if derivatives are 
absent) 



{(f>A{x),(l)B{y)} = lAB{x)5{x,y), 



(13) 



dF 



E°AimA. 



(8) 



This is fully justified if all variations d(l)A and all 
their spatial derivatives are zero on the bound- 
ary. In a more general case E'^(f) gives us only 
a part of the full variation. As a consequence 
of that fact Euler-Lagrange derivatives may not 



where 

lABix) ^ I^bDl, Ia^=I 



AB=^AB\4>^C^)- (14) 



The new feature of our formalism is that the sur- 
face contributions to this brackets are allowed Q 

{4>A{x),<t>B{v)} = 9iP{x)i^/J{x)S{x,y). (15) 
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For example, Ashtekar's transformation in the 
canonical gravity which is of the type (||) leads 
to the generalized form of the Poisson brackets 
given above [Q. This is a rather general fea- 
ture of transformations of this type. We con- 
sider another example connected with the non- 
linear Schrodinger equation in other place (l9| . 

Even if the Poisson brackets of fields do not 
contain surface contributions such contributions 
may arise in the calculations of the Poisson al- 
gebras for some transformation generators con- 
structed by means of these fields. It is so because 
these nonstandard terms may be a result of mov- 
ing the derivatives of the (5-function from one of 
its arguments to another. The standard rule 



iAB{x)S{x,y) ^ i\g{y)6{x,y), 



(16) 



is applicable but with the definition of the adjoint 
operator modified to preserve all the boundary 
terms 



= (-1) 



K\ (K -L 
M 



(17) 



X Dk-l-mIba 



For example, if we preserve the boundary contri- 
butions, then the usual formula 



should be replaced by the new one 



(18) 



enix) 



d 



0n{y)-^^]5{x,v)^-e^^h{x,y). 



This solves one paradox which arises in under- 
standing the result obtained for asymptotically 
flat spaces in the canonical General Relativity. It 
is shown (nJsl that 



(19) 



{i?(0,i?(r/)}«i?(K,ry]), 



where r]^ are the Killing vectors of the back- 
ground flat metric, rj\ is their Lie bracket and 
H{^), H(r]), H{[^, T]]) are generators with nonzero 
surface terms. The paradox is in the observation 
that if we first consider integrands of (n9[) and 



calculate their Poisson brackets according to (18) 



then we will get zero result due to the closed con- 
straint algebra of the General Relativity. 

The third step in our revision of the Poisson 
bracket formula is dictated by purely mathemati- 
cal reasons. Speaking in mathematical terms, the 
extension of the formal variational calculus pro- 
posed above is the introduction of a new grading. 

A grading in linear space L is a decomposition 
of it into a direct sum of subspaces, with a spe- 
cial value of some function p (grading function) 
assigned to all the elements of any subspace ||] . 

Here the function p takes its values in the set 
of all positive multi-indices and thus, 



L 



J=0 



(20) 



Elements of each subspace are called homoge- 
neous. 

A bilinear operation x,y i—> x o y^ defined on 
L, is said to be compatible with the grading if 
the product of any homogeneous elements is also 
homogeneous, and if 



p{xoy) ^p(x) +p{y). 



(21) 



It is necessary to define the pairing between 1- 
forms and 1-vectors, which then really will induce 
all other operations, as an operation compatible 
with the introduced grading. In our notions 1- 
forms are 



'n 



^AK^K 



whereas 1-vectors are 



Their bases are dual 
5 



Ha{x), 



5AB5{x,y), 



and the graded differential operators 
« = (^''nO^AlDK. 

4' = ^^[/Vi'ii^L, 



(22) 



(23) 



(24) 



(25) 
(26) 
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serve as coefficients of the decomposition over 
these bases. Let us define a trace as 

Tr(aV3) = e\^+'^DLaf^DK^^^l (27) 

that is evidently a biHnear and commutative op- 
eration, ft is also easy to check that 

TTiiDa)^)) = DTT{a{j}). (28) 

This important property of the trace operation 
allows us to use it for the definition of the pairing 

aiiP) ^i}Ja = J Tiiaip), (29) 

which is independent on the ambiguity in the rep- 
resentations of the operators a, V following from 
the freedom to do the formal integration by parts. 
For example, we can remove all the derivatives 
from the basis elements 

a = j e^^^afS^A, (30) 

* = (31) 

thus transforming a and -0 to the so-called canon- 
ical form (compare with ||^). This formal integra- 
tion by parts does not change any integral over 
the finite domain and is useful for the illustra- 
tion of analogy with the standard formalism. To 
return back to the usual formal variational calcu- 
lus we should only put 9n{x) = 1, i.e., Vl = M" 
then the "columns" like a^''^ will be reduced to 
their first terms a^*^^ and for the canonical repre- 
sentation of the 1-vector and 1-form their pairing 
will be reduced to the standard one 

V-Ja^l (32) 

After making the above three steps: the revi- 
sion of differential, Poisson bivector and pairing 
we can obtain the new formula. There are at least 
three ways to write it, in correspondence to the 
three ways to write the differential of a local func- 
tional: through the full variational derivatives 



through the Frechet derivatives 
{F,G} = j e[:^hT[f'J^^'^g's), (34) 
and through the higher Eulerian operators 
{F,G} = j e'^^Dp+Q{EP^{f)i^^lE%)). (35) 

3. EXAMPLES 

The calculation of the Poisson brackets by the 
new formulas can be made in not more compli- 
cated way than by the old ones. First, to get the 
Frechet derivative from the first variation is even 
easier than to get the standard Euler-Lagrange 
derivative because the integration by parts is not 
needed. Second, we can exploit covariance prop- 
erties and use the covariant derivatives instead of 
the ordinary ones. 

As an example, we calculate the Poisson brack- 
ets of the two spatial diffeomorphism generators 
in the canonical General Relativity for a finite 
domain 

H{N') = [ 7:'^VjN, + V,Nj)d^x^ (36) 
Jn 

= (f 2njN'dS.j - / 2N'Vjnfd^x. 
J do. Jn 

From the first variation we get 

/iU(7V'') = V,N, + V,N,, (37) 

and then calculate the Poisson bracket according 
to formula (|^ 

{H{N^),H{AP)}^H{[N,Mr), (38) 
where 

[N, Mf = N'^APj, - M''N\f,. (39) 

Let us mention that the calculation according to 
the standard formula gives additional surface con- 
tribution violating the diffeomorphism algebra 

AH = (f) ^T'^{VjNi + V^NJ)M''- 
Jan 

- {VjAU + V,Afj)7V'=) dSk. (40) 
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This violation is zero in the only case when 
and are the Killing vectors on the boundary. 
Therefore we have a free boundary closure of the 
spatial diffeomorphism algebra by means of the 
new brackets. 

As a second example let us consider the Yang- 
Mills field in a finite domain. It is suitable to use 
orthonormal curvilinear coordinates Xk in , so 
that they are compatible with the boundary d^l = 
{Xk : Xi = R = const, X — X2,X3}. If Xk are 
Cartesian coordinates then the local frame 

can be used. The Hamiltonian has the form ||2^ 

- A^,{d.E^-gt^''^A'^E^)^. (42) 

It should be accompanied by the surface contri- 
bution 

AHn = J^^ dX AS (J^^Ef,^ + x'^) , (43) 

where h = /ii/i2^3, and x"^ ^I'o the surface vari- 
ables possessing the Poisson bracket 

{X"(X), /(X')} = ge'' x'{±) S{X X'), (44) 

and commuting with the volume variables. After 
fixing the spatial Fock-Schwinger gauge, A°^^ = 
0, the Hamiltonian evolution on the boundary is 
given by 

E^,-,{R,±) = -<^^i?(\)(i?,X)AS(i?,X), (45) 
X^X) = <''=x''(X)Ag(i?,X), (46) 

where the boundary condition G"-^^ {R, X) = 0, 
Gf = '^jF^j compatible with the localized time 
evolution is also assumed. 

Here the Gauss law constraint is prolonged onto 
the boundary by introduction of the surface vari- 
ables x°(X). The standard approach requires 



a fixation of the static boundary conditions for 
E'^^j , and in its turn this requires that the La- 
grangian multiplier should be zero on the 
boundary. Then the boundary conditions are to 
be gauge-dependent or E'^^^ be zero. The ap- 
proach based on the dynamical boundary con- 
ditions permits to save gauge invariance on the 
boundary. 

The residual gauge invariance of the theory is 
manifested in the above dynamics on the bound- 
ary. This implies that the boundary conditions 
put onto E^j^s^ may not necessarily be arbitrary to 
preserve the gauge invariance. It has been argued 
that the dependence of the partition function on 
these boundary conditions may be considered as a 
confinement criterion in the SU(N) gauge theory 
| |2^ ] , and that the surface terms play an important 
role for understanding that phenomenon p^. 

4. CONCLUSION 

As we have seen from above, the Poisson struc- 
ture can be introduced prior to any boundary con- 
ditions. This is analogous to the Hamiltonian me- 
chanics where constraints are treated later than 
the Poisson brackets. We may expect that the 
treatment of the boundary conditions could pro- 
ceed similarly, so that primary and secondary, 
first and second class boundary conditions may 
arise. We may get an analog of the Dirac bracket 
at the end of the standard reduction procedure. 

What are the Hamiltonian equations generated 
by the new bracket? They can be called as a weak 
form of the equations of motion If we try 
to understand ^Q-functions as distributions seri- 
ously, then we will have singular boundary terms 
in the equations and may encounter with ambi- 
guities in solving such equations. Therefore, the 
construction of the closed Poisson algebras with 
surface terms seems to be a more promising direc- 
tion. There we deal with local functionals, rather 
than functions, and the Poisson bracket does not 
move us out of that class. It is quite possible 
that the Hamiltonian dynamics for the function- 
als may become of more importance for the quan- 
tum field theory than for the classical one. 
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